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ABSTRACT 

The amount and distribution of heavy elements in Jupiter gives indications on the process of its 
formation and evolution. Core mass and metallicity predictions however depend on the equations of 
state used, and on model assumptions. We present an improved ab initio hydrogen equation of state, 
H-RE0S.2 and compute the internal structure and thermal evolution of Jupiter within the standard 
three-layer approach. The advance over our previous Jupiter models with H-REOS.l by Nettelmann 
et al. (2008) is that the new models are also consistent with the observed > 2 times solar heavy element 
abundances in Jupiter's atmosphere. Such models have a rock core mass Ale = 0-8 M®, total mass 
of heavy elements Mz = 28-32 M^, a deep internal layer boundary at > 4 Mbar, and a cooling time 
of 4.4-5.0 Gyrs when assuming homogeneous evolution. We also calculate two-layer models in the 
manner of Militzer et al. (2008) and find a comparable large core of 16-21 M^, out of which ^ 11 
is helium, but a significantly higher envelope metallicity of 4.5 x solar. According to our preferred 
three-layer models, neither the characteristic frequency {vq ^ 156 /iHz) nor the normalized moment 
of inertia (A-^ 0.276) are sensitive to the core mass but accurate measurements could well help to rule 
out some classes of models. 

Subject headings: planets and satellites: individual(Jupiter) - equation of state 



1. INTRODUCTION 

Numerous discoveries of unusual hot Jupiters, in 
particular those that appear to be abn ormally in- 
flated (HD 209458b; 
metal-rich (HAT-P-2b 
rupted (WASP-12b; E 



Baraffe et al.l 



2005 ) , extremely 



Leconte et al.l[2009t ). tidally dis- 
et all l2O10[r that are thought 



to have formed by grav itational instability (HR 8799b- 
e; iBaruteau et al.l 120 lit) , or that even allow to study 
non-equilibriu m chemistry in ir radiated atmospheres 
(HD 189733b: lR)rtnev et al.ll2010[ ). have drawn some at- 
tention away from Jupiter. 

While models for exoplanets often span a variety of 
possible solutions due to the small number of observ- 
ables characterizing the planet such as radius, mass, age, 
and their rather large error bars, with Jupiter we face 
the opposite problem that many models fail to repro- 
duc e the available observati onal constraints (for review, 
see ISaumon fc GuillotI 1200^ . With the upcoming Juno 
mission, Jupiter will not only be back in focus, but the 
new data may also invalidate models that cannot cope 
with the expected accurate measurements for the grav- 
ity field, the moment of inertia, or the oxygen abundance 
in the atmosphere. 

As a major source of influen ce on Jupiter mod- 
els, iGudkova fc Zharkovl (Il999b0 identified the dis- 
tribut ion of metals an d heli um onto certain lay- 
ers. ISaumon fc Guillo^ (I2004D demonstrated in ad- 
dition the sensitivity of the Jupiter models to the 
equation of state (EOS) of hydrogen. In the past 
few years, the H EOS in the warm dense matter 
region where hydrogen metallizes has been improved 



ICaillabet et al.|[2"011[ ). The corresponding Jupiter mod- 
els were consistent with the avail able observational con- 
straints (jNettelmann et al.l [20081 hereafter N08). How- 
ever, the measured enrichment of Jupiter's atmosphere 
with noble gases (~ 2x solar), and with carbon and ni- 
trogen (~ 4x solar) suggests that oxygen is equally abun- 
dant, leading to an atmosphere metallicity of 2-4 x solar. 
Such high values could not be obtained with the for- 
mer H EOS (H-REOS.l) used in N08, which instead re- 
quired solar metallicity atmospheres. This low atmo- 
spheric metallicity was confirmed by another group that 
used the same met hod of com puting the high-pressure 
EOS data (Militze r et al.l [200^ . Apart from this sim- 
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ilarity however, the Jupiter models suggested by these 
two groups differed from each other substantially with 
respect to core mass, heavy element content, and inter- 
nal layering. In this paper we present new Jupiter mod- 
els with ^ 2x solar metallicity atmospheres as a result 
of applying an improved ab initio H EOS. 

In §[l]we compare the improved H EOS (H-RE0S.2) 
with H-REOS.l and describe our procedure of structure 
and evolution modeling. In §[3] we present the core mass 
and metallicity of our new Jupiter models and show that 
the d ifference in the Jupiter models of IMilitzer et al.l 
(j2008f ) and our group does not originate from the as- 
sumption of two or three layers (§ 13. 2p . We also present 
calculations of the moment of inertia (5 13. 3|) . of the char- 
acteristic frequency i^o of the global oscillations fi) l3.4p . 
and of the thermal evolution (§113]). We discuss Jupiter's 
atmospheric metallicity in § |4] and give our conclusions 
in §[5l Appended to this paper are descriptions of our 
method to calculate the entropy (Appendix |A| and the 
moment of inertia (Appendix IB|) . 

2. METHODS 
2.1. Improved H EOS 
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Our new hydrogen equation of state (H-RE0S.2) is 
constructed in the same way as the former version H- 
REOS.l (N08). The H EOS is assembled from different 
contributions at the low density range [p < 0.1 gcm~'^) 
and the high density range (0.2 < p < 9 gcm-3). At 
low densities, pressure dissociation and pressure ioniza- 
tion do not play an important role for the EOS. For these 
densities, models within the chemical picture agree well 
with experimental data. As for H-REOS.l, w e here apply 
the F luid Variational Theory (FVT+) model f Hoi st et al.] 
I2007D where ionization is neglected for temperatures be- 
low 3000 K. 

In the high density range, our EOS table is ob- 
tained from finite-temperature density-functional the- 
ory molecular dynamic s (FT-DFT-MD) simulations 
using the code VASP (iKresse fc Hafned [l99a . .1994; 
iKresse fc Furthmiilled 11996( 1. There are three improve- 
ments of our new H EOS. First, due to the enormous 
increase in computing power we were able to perform 
the simulatio ns with 256 partic les in a box compared to 
64 particles (iHolst et al.l 120081 ). Second, we performed 
them up to 50,000 K, in order to cover the interior of a 
warm, young Jupiter, compared to 20,000 K (jHolst et al.l 
120081 N08). In the simulations, a quantum mechanical 
treatment of the molecular vibrations of the ions is not 
included. Thus thirdly, after completion of the simula- 
tions, we add the energy of a quantum mechanical oscil- 
lator, Umb = fcB6m&(0.5 -I- {exp{Qmb/T) - 1)"^), to the 
internal energy of each molecule and subtract its classical 
energy of vibration, k^T. The degree of di ssociation is 
estim ated via the coordination number (see lHolst et al.l 
|2008[) . This contribution is non-negligible if the temper- 
ature o£_thfi_s;Y^tem_2S_bel^ the vibration temperature 
&y,b (|French fc Redmedl2009( ) . For hydro gen, we use the 
experimental value &vib = 6338.2 K. 

For H-REOS.l we obtained convergence of the EOS 
data within 3% and estimated the overall uncertainty 
due to both statistical and systematic errors to be below 
5%. The new ab initio data in H-RE0S.2 have an uncer- 
tainty of 1% for most of the density-temperature space. 
They are within the uncertainty of the former data for 
p > 0.3 gcm"-^ but shifted to 2% higher pressures in a 
T-p region where pressure dissociation and ionization oc- 
curs. Figure [T] shows selected isotherms. With the onset 
of pressure dissociation at 0.2 gcm"'^, the difference rises 
to 10%. This is a rather small density for FT-DFT-MD 
simulations that requires large computational effort to 
reach converged results. The displayed region of 5000- 
10000 K and 0.2-1 gcm~'^ is relevant for t he Jupiter adi- 
abat. Here, the SCvH-i hydrogen EOS (jSaumon et al.l 
|1995[ ) shows systematically higher pressures. At densi- 
ties of 0.1 < p < 0.2 gcm"-^, the difference in pressure is 
a consequence of a revised interpolation between the low 
density regime and the data points at 0.2 gcm~'^. This 
difference is particularly pronounced for the internal en- 
ergy (not shown). Along the Jupiter adiabat (Fig. [2]), 
the revised internal energy requires higher temperatures 
in order to keep the entropy constant at the level defined 
by Jupiter's outer boundary condition of 170 K at 1 bar 
(see ij |2.3p . and also leads to a decrease in the density at 
low pressures of 0.1 to 10 GPa. However, due to little 
mass prevalent at those pressures in Jupiter, this region 
does not affect our Jupiter interior models. Instead, it is 



the 10-50 GPa pressure region where the lower densities 
of H-RE0S.2 map onto the new Jupiter models. 

Our EOS data tables provide the thermal EOS P{p, T) 
and the caloric EOS u{p, T), but not the specific entropy 
s(/9, T). In Appendix [Al we describe our method to derive 
an adiabat through a given point (P, T) . 

2.2. Code improvements 

Careful inspection of our code used for the calculation 
of Jupiter's shape and rotation in N08 revealed that two 
coefficients, (one of second-order, the other one of third- 
order) contained errors. Their effect on Jupiter models 
with H-REOS.l is twofold; i) the layer boundary shifts 
from ~ 8 to ^ 4 Mbar when the outer envelope is made 
to have solar metallicity; ii) the metallicity increases by 
~ 0.02 (in absolute value) in the outer envelope and 
~ 0.035 in the inner one so that the maximum total mass 
of metals increases from 32 to 41 M^, where the latter 
value was presented in N08. All other results are almost 
unchanged. The two errors were found by a comparison 
of the resulting gravitational coefficients J2„, n = 1-4, of 
linear and polytropic density Jup iter and Saturn mod- 
els wit h analytic solutions given in IZharkov fc Trubitsvnl 
()1978() . and with numerical solutions for these density 
distributi ons by T. G u illot ( pers. comm. 2008). For de- 
tails, see iNettelmannI ()2009[ ). Another improvement of 
the code is the application of a Newton-Raphson scheme 
to find the central conditions that meet the envelope 
pressure and mass at the core-mantle boundary. For the 
current version of our code, we here introduce the la- 
bel MOGROP-11 (MOdellierungsprogramm fiir GROsse 
Planeten, 2011.) 

2.3. Interior structure and evolution modeling 

Apart from the code improvements mentioned above, 
our Jupiter interior and thermal evolution models are 
generat ed by essent i ally t he sa me code as introduced 
in N0 8, INettelmannI (|2011l ). and iFortnev fc NettelmannI 
(|2010( ). where in the latter case the code was applied to 
the structure and evolution calculations of Uranus and 
Neptune. We here recall the properties of our three- 
layer structure model for giant planets, review the ob- 
servational constraints for the new Jupiter models, and 
describe the method of structure and evolution modeling. 



Three-layer structure model. — A sufficient number of 
parameters to meet the observational constraints and 
boundary conditions are available if we assume a three- 
layer structure with two envelopes and a core. The en- 
velopes are adiabatic, homogeneous, and consist of hy- 
drogen, helium, and heavy elements, while the core is 
made of rocks. Heavy elements in the envelopes are 
represented by the water equation of state H2O-REOS 
' French et al]"l20ql N0 8'). for helium we use He-REOS 



..Kietzmann et al.l2007l: N08 ) , and for rocks the formula 
given m IHubbard fc MarlevI (|1989f ). The outer envelope 
has a lower He abundance (Yi) than the inner envelope 
(12)- We allow for different heavy element mass fractions 
in the outer {Zi) and in the inner envelope (.^2)- The 
transition pressure between the envelopes Pi_2 is a free 
parameter. The density and entropy are discontinuous 
at layer boundaries. 
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Observational constraints. — The models are required to 
meet the observational constraints for the total mass Mj , 
equatorial radius at the 1-bar pressure level 1-bar 
temperature Ti, solid-body period of rotation w, grav- 
itational moments J2, J4, Jei atmospheric helium abun- 
dance Yatm, and have the same He mass fractioiE Y as 
the protosolar cloud once had. As long as the oxygen 
abundance below Jupiter's water cloud deck is not mea- 
sured, only the lower limit for the heavy element abun- 
dance Ziitm in Jupiter's atmosphere is known. This value 
is 1 X Zq (see jj IH) , where Zp, is the protosolar cloud's 
metallicity of 1.5% (|Loddersll200"l . 

These observational constraints for the structure mod- 
els can be divided into four groups depending on whether 
they enter the modeling procedure explicitly or need to 
be fitted by adjusting model parameters, and whether 
or not they are varied within the la error bars. The 
first set (explicit, varied) is empty. The second set 
(explicit, not varied) consists of Mj = 317.8338 M^, 
2tt/lu =9h 55m 30s (the period of the magnetic field), 
the upper limit Ti = 170 K, and Fatm = 0.238. While 
the small uncertainties in Mj and uj would not affect the 
resulting models, lowering Ti by 5 K would reduce Zi 
by ^ 0.3 Zq (N08), in the opposite direction to what we 
are aiming for. Modifying Yatm by Icr (±0.007) changes 
Zi by about =f0.3 Zq. Because Yatm was measured by 
the Galileo entry probe at a depth where vertical convec- 
tion ensures homogeneous distribution of those elements 
that do not form clouds (such as the noble gases) deeper 
inside, we set Yi = Yatm- The parameters of the third 
group (to be fitted, not varied) are Rj = 7.1492 x lO^m 
and Y = 0.275, whereof i?j is fitted by choosing the mean 
radius i?,n of the equipotential surface that coincides with 
the 1-bar pressure level, and Y by adjusting Y2. We find 
i?,„ = 10.955 i?e and Y2 = 0.285-0.325. The gravita- 
tional moments J2, J4, and Jg build the fourth group (to 
be fitted, varied). In particular, Z2 is used to reproduce 
J2/10~^ within the tight observational bounds 14697(1), 
and Zi is used to adjust J4/10~^ to either —587 (the ob- 
served central v alue) or —589 (the observed upper limit) 
()Jacobsonll2003[ ). When J2 and J4 are matched, Je/lO"^ 
is found to lie within 34.0 and 37.5, consist ent with the 
observed limits (29.1-39.5; I Jacobsoiil 120031 ). Our com- 
putation of Jupiter's thermal evolution makes use of the 
observed effective temperature T^s = 124.4 ± 0.3 K. 

Modeling procedures. — Internal profiles of m, P, T, and 
p which depend on the radial coordinate I are obtained 
by numerically integrating the structure equations 

p~^dP/dl = -d{V + Q)/dl (1) 

and 

dm/dl =4:nl^ p . (2) 

Equation ([T]) states that the force acting on a mass ele- 
ment with density p at radius / is balanced by a pressure 
gradient, where the force —d{V + Q)/dl arises from the 
gravity field V and the centrifugal force with potential 
Q due to rigid-body rotation. Equation ([2]) expresses 
the amount of mass dm included in a spherical shell of 

^ Throughout this paper, the label Y denotes a He mass fraction 
with respect to the H-He subsystem. 



density p and thickness dl at radius I. These two ordi- 
nary, first-order differential equations are integrated in- 
ward starting at the surface / = with the two bound- 
ary conditions P(i?m) = 1 bar and m{R-^) = In 
order to satisfy in addition the inner boundary condi- 
tion to(0) = 0, a further parameter is needed. For this 
we choose the core mass M^. Thus for given values of 
the observational constraints and of the transition pres- 
sure Pi_2, the result is one single Jupiter model, ob- 
tained in terms of the internal profile and values for Zi, 
Z2, and Mc. To calculate the planetary shape and the 
gravitational moments we apply the Theory of figures 
(jZharkov fc Trubitsvnlll978D up to third order as in N08 
and also to fourth order. Such structures describe Jupiter 
at present time io = 4.56 Gyr. 

For modeling Jupiter's t hermal evolution we apply a 
procedure similar to that of lSaumon et al.l (jl992l) . To de- 
scribe Jupiter's interior at earlier times t Kt^, we select 
a few representative structure models and generate for 
each of them a sequence of ~ 80 warmer interior profiles 
by increasing Ti up to 800 K while m(Pi_2), Mc, and the 
clement abundances are kept at constant values. Further 
models with intermediate Ti values arc then generated by 
interpolation. While for the structure models at present 
time with known R]{to) the core mass was chosen to en- 
sure mass conservation, we here invert the problem to 
find the planet radius R m H) for given core mass. Un- 
like [Saumo^^F^l] (jl992l) . rotation is included in the ap- 
proximation of spherical symmetry where the centrifugal 
force reduces to the zero-order term —dQ/dl = (2/3)a;^/. 
We calculate the cooling curves in the usual approxi- 
mation of constant angula r velocity (jSaumon fc Guillo^ 
120041: iFortnev et al.|[201l[) . and then also by in cluding 
angular momentum conservation (|Hubbarcill 19701 ) , which 
requires just 1-3 additional iterations for a given value 
of Ti , and the corresponding change in the energy of ro- 
tation, dEj-of The time dt passed between profiles with 
different surface temperatures and hence different inter- 
nal entropies is given by the equation of energy balance, 

LcS — Lcq — —dEint/dt , (3) 

where Leff is the planet's observable luminosity in the 
infrared and Loq the luminosity the planet would have 
in equilibrium with the insolation. Using the Stefan- 
Boltzmann law for the energy emitted by a black body 
with uniform temperature, these luminosities can be used 
to define the commonly used temperatures Teg, Toq, and 
Tint. T his way. Teg (t p ) is derived from the Leq value 
given in iGuillotI (|2005[ ). It is set constant over time or, 
alternatively, varied with time according to a linear in- 
crease of the irradiation flux := Leq/47ri?^, start- 
ing with Fcq(O) = 0.7Foq(to) as predicted by theoretical 
standard models for the Sun (IBahcall fc Pinsonneaultl 
[19951 ISackmann fc Boothrovdl I2003D . Note that the 
solar standard evolution model predicts an increasing 
bolometric luminosity with time, whereas the short- 
wavelengths (chromospheric and coronal) emissions of 
the young Sun may have been up to 1000 times stronger 
than those of the present Sun as indicated by the mea- 
sured XUV fluxes of young solar-like stars (iRibas et al.l 
I2005D . Absorbed XUV irradiation can lead to heating 
of and mass loss fro m the upper planetary atmosphere 
(jLammer et al.|[2003[ ) but absorption at high altitudes is 
not suspected to influence the energy balance of the in- 
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terior (|Guillot fc ShowmaD|[200^ . Therefore, we Ignore 
the activity-related, time-dependent solar XUV flux in 
our evolution calculations for Jupiter. 

The difference between the luminosities radiated away 
and absorbed equals the loss of the planet's intrin- 
sic energy dE'int per time, see Eq. ([3]). In the enve- 
lope, the energy lost per mass she ll simp ly is the heat 
6q = T{m) ds{m). IGuillot et all (fT995l) offer a con- 
venient closure relation by relating T^s to Ti, namely 
Ti = K Tli ^^ q-0.167 acco rding to the model atmo- 
sphere grid bv lGraboske et al . (1975), where the param- 
eter K can be used to reproduce the obs erved Tpff(tn), 
and q is the surface gravity. Also like IGuillot et al.l 
we take into account the core's contribution to 
the planet's intrinsic luminosity. The heat loss of the 
core is c,, rfTf^nrs , where we use Ci, = 1 J g"-'^ as in 
IGuillot et al.l (jl995| ) , and dTcorc is the temperature differ- 
ence of the isothermal cores of subsequent interior profiles 
with different T\ values, as is ds{m) the entropy differ- 
ence at mass shell m between those interior profiles. Al- 
though the abundances of radioactive elements decrease 
exponentially with time so that the energy production 
from radioactive decay in rock with meteoric composi- 
tion may have been an order of magnitude larger in the 
past than it is today, we here use the present Earth's 
value iradio = 2 X 10^'^ Js^^M^^ to account for this 
contribution. This treatment seems justified for Jupiter 
because we find the prolongation of the cooling time due 
to the core's heat loss to be O.OI Gyr only. Thus the 
cooling equation reads 

J 1^1^' dmTds + McCydTcoie + dE^-ot 

After integrating Eq. (j4]) backward in time we obtain the 
cooling time r that the particular Jupiter model needs 
to cool down from an arbitrarily hot initial state to the 
present state. 

3. RESULTS 

3.1. Core mass and metallicity of three-layer models 

Results for the core mass, the outer envelope metal- 
licity, and the inner envelope metallicity as functions of 
the depth of the layer boundary between the envelopes 
are shown in Fig. El If the layer boundary is located 
higher in the planet (lower transition pressures), Zi and 
Z2 become smaller. The lighter envelopes then require 
a larger core mass to conserve the total mass. When 
Zi decreases down to zero, the core mass cannot grow 
any further. However, we do not consider models with 
Zi < Zq acceptable as this is inconsistent with Jupiter's 
observed atmospheric abundances. We find Z2 S> Zi 
and Mc < 10 for all our acceptable three-layer mod- 
els. As the gravitational moments are most sensitive to 
the density in the outer part of the planet at pressures 
of a few Mbars ( J2) or even below 1 Mbar (J4) and in- 
crease with the local density, Zi must decrease when the 
denser inner envelope extends farther out. Otherwise, 
I J4I would become too large. On the other hand, as we 
go to higher pressures deeper inside the planet, the sensi- 
tivity of J4, and later also of J2 approaches zero. There- 
fore, Zi changes only weakly with P1-2 for deep internal 
layer boundaries. However, since J2 is adjusted by Z2, 



Z2 must rise strongly with P1-2 in order to provide a suf- 
ficiently high mass density that is able to keep J2 on the 
high level of the observed value. With rising metallicity 
the envelope becomes denser, so the core mass must de- 
crease. When Mc = is reached, the layer boundary is 
as deep as it can be, and Zi adopts a maximum value. 
For models with H-REOS.l, this maximum is 1.0 Zq if 
J4/10-^ is adjusted to 587, and 1.5 Zq if 74/10"^ is ad- 
justed to 589. For models with our improved H EOS, 
this maximum is 2.5 Zq (J4/10"^ = 587) and 2.7 Zq 
(J4/10^^ = 589), respectively. This enhancement in Zi 
arises from the ^2% higher pressures at given hydro- 
gen mass density (Fig.[T]) , implying ~ 2% lower hydrogen 
mass density at given pressure. The lower partial density 
of hydrogen in the mixture requires to be compensated 
for by a similar amount of metals. 

Including 4th order coefficients to the calculations of 
Jupiter's shape and rotation systematically shifts the so- 
lutions to ^ 0.8 Zq lower outer envelope metallicities, 
inducing an additional uncertainty on Zi that is of the 
same size as the error bars of single observables such as 
Ti and J4, see below. For Mc and Z2, the inclusion of 
higher-order terms has an even bigger effect than the 
improvement in the H EOS. Obviously, a convergence 
check of the solution with increasing order of the theory 
is highly important but beyond the scope of this paper. 
The accuracy of our current computer code is insufficient 
to reliably study the effect of higher than 4th order co- 
efficients. From preliminary calculations, solutions with 
5th order coefficients are found to lie between those of 
3rd and 4th order. 

Figure [3] also shows the observed element abundances 
in Jupiter's atmosphere. Assuming that O is equally en- 
riched as C, these abundances indicate an average enrich- 
ment of 2 (noble gases) to 4 (C,N,0). Taking 2x solar 
as the lower limit for Zi, models with P1-2 < 4 Mbar 
drop out of the realm of acceptable models, reducing the 
upper limit for the core mass down to 8 Mq. The total 
mass of metals is Mz = 29-32 Mq. 

As representatives for our acceptable models we rec- 
ommend models Jll-4a and Jll-8a, highlighted in Fig.|31 
Model Jll-4a is the one with the outermost possible layer 
boundary, i.e. at 4 Mbar. Model Jll-8a has P1-2 ~ 
8 Mbar for which Zi is > 2Zq for all of the above 
considered uncertainties. A machine-readable table for 
model Jll-4a is provided as supplemental material, see 
Table [TJ For these two models we also vary Y, Yi, Ti, 
and J4 within their 1-cr error bars. Modifying Y by 
Ay = ±0.05 influences I2 in the same direction, and 
thus leads to a change in Z2 of ^8%. Modifying Yi by 
AYi = ±0.05 causes a response in Zi of =f10%. Similar- 
ily, a 5 K colder 1-bar level leads to a 50% decrease in 
Zi. In both cases, Z2 consequently changes by ±3-4% in 
the opposite direction to Zi, and Mc changes by =f2-3%, 
again opposite to Z2. The colder outer envelope adiabat 
also cools the interior, resulting in a core that is 100 K 
colder. The slightly denser H-Hc subsystem that results 
reduces the amount of metals in the envelope required 
to match J2, but the net effect of ATi on Z2 remains 
positive. Setting |J4|/10~^ down to 585 lowers Zi by 
^ 15%. Clearly, the response of Z2 (increase) and thus 
of Mc (decrease) in this case is stronger than in case of 
ATi and Ali. Including these uncertainties, the total 
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mass of metals changes but insignificantly to 28-32 M^. 
3.2. Core mass and metallicity of two-layer models 
The striking difference between the iMilitzer et al.l 



(j2008[ ) (hereafter M08) predictions for the core mass of 
Jupiter (16 ± 2 M®) and that of our group (0-8 M©), 
where both groups applied similar ab initio simulations 
to generate the EOS data, has pr ompted speculations 
abou t causes for that difference (jMilitzer &: HubbardI 
|2009| ). We here investigate how this difference relates to 
the modeling assumptions of two-layer (M08) or three- 
layer (N08, this work) structures. For that purpose we 
have calculated two-layer models following M08. As in 
M08, the homogeneous envelope has Yi ~ 0.238, and the 
missing helium to obtain an average Y of 0.275 is not ex- 
plicitly accounted for but may be part of the core mass. 
The core is either pure rocks or rock-ice, where we switch 
from ice to rock at 70 Mbar. As a result, only the inner- 
most 4 Mq are rocks. Because these models do not have 
the degree of freedom Z2, we can only match either J2 
or J4. As in M08 we choose to match J2 by varying the 
envelope meta llicity Z. We apply our usual numerical 
procedure (see iNettelma nn 2011 for details) and vary Z, 
starting from 0.01, until J2{Z) meets the observed value. 
The core mass must decrease with increasing envelope 
metallicity (denser envelopes) in order to ensure total 
mass conservation. If the number of layers assumed for 
the structure of the model were responsible for the differ- 
ent Jupiter models we would expect to obtain the same 
two-layer model as in M08. This is not the case, as illus- 
trated in Fig.H While the M08 models require 4 ± 2 
heavy elements in the envelope to meet J2, correspond- 
ing to Z = 1.0 ± 0.5 Zq for a 16 Mq core, our two- layer 
models require Z = 4.3 Zq. 

This significantly larger envelope metallicity of our 
two-layer models enables us to make a transition to 
three-layer models where metals are shuffled from the 
outer part of the envelope to the inner part in order get 
I J4|/10~^ down to the observed value (from 611 to 587). 
The M08 models do not have this degree of freedom as 
the already low envelope metallicity would become zero 
long before J4 is matched. Hence the M08 models stay at 
J4/10~^ = 614 and a large core while our two-layer mod- 
els can undergo a transition to three-layer models where 
Zi becomes smaller and Z2 larger than in the homoge- 
neous envelope case. As a consequence of large Z2 values, 
the core mass decreases further. In contrast, the core 
mass of the two-layer models is large because this mass 
also contains the missing mass of helium dMno = (0.275- 
0.238)(A/j " Mc) ~ 11 Me needed for an overah abun- 
dance Y = 0.275. Thus the ice-rock core mass of two- 
layer models is about 5-10 Mq, just at the upper limit 
of our three-layer models. This estimate would also hold 
for the two-layer model of M08. However, we cannot 
reproduce their Jupiter model, the reason for which re- 
mains unexplained at the moment. We conclude that 
the difference between the M08 two-layer models and our 
three-layer models does not originate from assumptions 
about the number of layers, but instead is a symptom of 
already different models within the two-layer frame. 

3.3. Moment of inertia 

We have calculated the moment of inertia / as de- 
scribed in Appendix |BJ The non-dimensional form (A) 



is obtained as A = I Rm^- For all of our three- layer 
models, we find an almost invariant value A = 0.27605 ± 
0.03%. If we ignore Jupiter's shape deformation, the 
resulting A decreases by 4% down to 0.26539 ± 0.03%. 
This is clos e to t he predictions of 0.2629—0.2645 of 
iHelled et al.l (j2011l ). who considered a broad set of in- 
terior models with core masses between and 40 Mq as 
allowed by the use of six-order polynomials to represent 
the pressure-density relation in Jupiter's envelope. From 
our calculations we conclude that A is not an appropri- 
ate parame ter to constrain Jup iter's core mass further. 
In contrast. IHelled et al.l (|2011| ) found a significant vari- 
ation of 0.6% over the full range of their models and still 
a 0.2% variation when excluding models with Mc larger 
than 10 Mq. Juno's measurement of Jupiter's A is ex - 
pected to have an accuracy of 0.2% (jHelled et al.ll20l"ll ). 
hence sufficient to distinguish between the different pre- 
dictions. 

3.4. Characteristic frequency 

Jovian seismology has been long recognized as a 
unique opportunity to infer interior structure properties 
(jVorontsov et al.iri976( ). In particular, the predicted low- 
frequency acoustic free oscillations depend sensitively on 
the core mass and i nternal layer boundaries of t heoreti- 
cal Jupiter models (jGudkova &: Zharkovlll99^a|) . How- 
ever, the expected small amplitudes (~ 10 cm) and 
Jupiter's rapid rotation m ake it difficult to m easure 
acoustic modes. Recently, iGaulme et al.l (|2011| ) used 
the SYMPA spectrometer to detect the lowest frequency 
modes 1000 fiRz) of Jupiter from spatially resolved 
radial velocity measurements. In the approximation of 
low degrees I and overtones of high radial order n ^ I, 
the acoustic modes f„_; are proportional to the charac- 
teristic frequency vq (also called equidistance) and are 
equally spaced by for given /: Vni ^ [n + 1/2)vq. 
IGaulme et all (pOlTI ) observed = 155.3 ± 2.2 /iHz. The 
inverse of vq is twice the time a sound wave needs to 
travel from the center to the surface (Jupiter's tropo- 
sphere). 



1^0 



2 / drc-^ 



(5) 



a useful parameter that interior models can easily be 
compa red to. In Eq. [5l c.^ is the adiabatic sound ve- 
locity. iGudkova fc Zharkovl (jl999al ) calculate vq = 152- 
155 /iHz for SCvH-i EOS based Ju piter models with five 
layers and core masses of 3-10 M^. iGaulme et al.l ([201TI) 
obtain the same range of values for vq with three-layer 
models of Jupiter that have a small internal He discon- 
tinuity and core masses of 0-6 . For our three-layer 
models with Mc = 0-8 M0 we find i'q = 155. 7-156. 3/.iHz. 
Since this range is rather narrow and slightly above that 
of former Jupiter models, a more accurate determination 
of the global mode spacing may help to rule out some 
Jupiter models but not constrain Jupiter's core mass fur- 
ther within the current uncertainty of ~ 0-10 M^. 

3.5. Cooling curves 

For a better comparison with published calculations of 
the cooling of Jupiter we first assume a constant irradi- 
ation fiux and constant angular velocity over time. The 
latter assumption is justified by the rather small change 
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of the rotational energy durin g evolution com pared to 
Jupiter's intrinsic energy loss (jHubbardl I1977D and the 
small change in planet radius over m ost of the contrac- 
tion time scale (jGuillot et alj|1995f ). Our H-RE0S.2 
based Jupiter model Jll-4a has a calculated cooling time 
of 4.66 ± 0.04 Gyr, see Fig. El and model Jll-8a of 
4.68±0.04 Gyr, where the uncertainty mostly arises from 
the observational error bar of T^ff. This result can be con- 
sidered in good agreement with the age of the solar sys- 
tem (tq = 4.56 Gyrs) given that none of the alternative 
proposed Jupiter models that rely on free-energy mod - 
els for the equation of state (jSaumon fc Guillotl l200l . 
such as LM-H4 (4.0 Gyr), SCvH-i (4.7 Gyr), and LM- 
SOC P (4.8 Gyr) gives a be tter agreement. On the other 
hand, iFortnev et al.l ()2011l ) investig ated the unce rtainty 
arising from the use of the iGrabosk c et al.l (jl975| ) model 
atmosphere grid for Jupiter's evolution. If instead a self- 
consistent model atmosphere grid was applied as it is the 
common approach for exoplanets, Jupiter's cooling time 
increases by ^ 0.5 Gyr. Sinking He droplets from H-He 
phase separation in Jupiter would also prolong the cool- 
ing time further. A process working in the opposite direc- 
tion could be rising material from core erosion. During 
the first 500 Myr, Jupiter has shrunken from 1.4-1.5i?j 
down to (Fig. [5]). Because the initial conditions 

for the long-term e volution of a 1 Mj planet are forgot- 
ten after 0.01 Gyr (iMarlev et al.l 120071 ) . we consider the 
displayed radius evolution at young ages realistic. 

In a second step, we keep the angular momentum con- 
served to Jupiter's current value (within a numerical 
accuracy of 0.4% which gives sufficiently smooth evolu- 
tion curves) and include the subseque nt change of rota - 
tional energy dE^ot{t), as also done in iHubbardI ()1970[ ). 
Such second order effects are important for estimates of 
the size of the correction factors that are necessary to 
bring the cooling time in agreement with tq. It is clear 
that both effects (A and diSrot) will prolong the cool- 
ing time: i) for the same total mass, a larger planet 
(the young Jupiter) rotates slower, implying a weaker 
centrifugal force that pushes matter outward, and the 
consequently smaller radius will allow less energy to be 
radiated away from the surface; ii) angular momentum 
conservation (L ~ luj) then implies a lower energy of ro- 
tation i^rot = 1/2 Luj at young ages. The increase in E^ot 
with time must be compensated for by a reduced lumi- 
nosity, implying again a longer cooling time. For model 
Jll-4a, the first effect increases r by 0.1 Gyr, and the 
second one by 0.2 Gyr, so that we find t = 4.96 Gyr. In 
a third step, we consider a time-dependent solar irradi- 
ation by a linear approximation of the Sun's luminosity 
evolution, which is assumed to have started with 70% 
of the current value. We find the lower irradiation to 
speed up the cooling by ^ 0.6 Gyr, so that we end up 
with r = 4.41 ± 0.04 Gyr. Figure [6] shows the relations 
between the structure parameter Ti that determines the 
planet radius, P = 27r/w and A. Their evolution with 
time can be read from the cooling curve Ti{t). 

The resulting cooling times of 4.4 Gyr (linearly in- 
creasing insolation) to 5.0 Gyr (constant insolation) sug- 
gest a reasonable understanding of Jupiter's interior and 
evolution, where remaining uncertainties seem to be at- 
tributable to uncertainties in the luminosity evolution of 
the Sun. This may be counter-intuitive as the present 
Sun is, by means of accurate helioseismology, far better 



constrained than Jupiter. For instance, uncertainties in 
the sound velocity profile of t he solar standard model 
are o f the order of 0.1% only (jBoothrovd fc SackmannI 
I2003D . On the other hand, non-standard solar evolution 
mode ls that predict a b right and more massive young 
Sun (jSackmann k, Boothrovd .20031 ) cannot completely 
be ruled out by observ ationally derived stellar mass loss 
rates, see lGiidell (|2007[) for a review. 

However, homogeneous, adiabatic evolution models for 
Saturn and Uranus that are equivalent to those presented 
here for Ju piter keep failing to r eproduce the observed lu- 
minosities (IFortnev et al.|[20TTh . indicating that the stan- 
dard three-layer model assumption may be too simplistic 
for so me giant planets. Therefore, iLeconte fc Chabrierl 
()2012[) investigated the possibility of layered convection 
in Jupiter, where heat is transported inefficiently by dif- 
fusion across thin stable layers, separated vertically by 
convective cells. They determined the maximum su- 
per adiabaticity that would give density distributions 
in agreement with the gravity field data. The deduced 
presence of roughly 10^ diffusive interfaces in Jupiter will 
qualitatively reduce the heat fi ux out o f the deep interior 
and shorten the cooling time (jStevenso n 198^. Quan- 
titative estimates of (at least) this effect are necessary 
before one can conclude a reasonable understanding of 
Jupiter's evolution. 

4. DISCUSSION 

According to our new Jupiter models with the im- 
proved H EOS (H-RE0S.2), metals are enriched by a 
factor of at most 2.7 in the outer envelope and atmo- 
sphere. That value can be increased to 3.0 if the atmo- 
spheric He abundance is lowered by Icr, but remains close 
to the lower limit of the measured C and N abundances. 
However, Jupiter's atmospheric metallicity is not directly 
observable. What can be measured are the abundances of 
single chemical species from which the atomic ratios with 
respect to the number of H atoms can then be inferred. 
Such measurements have been obtained for Jupiter by re- 
mote sensing with the Voyager IR and UV spectrometers, 
IR spectrometry with the Earth Orbiter ISO, Galileo or- 
biter measurements in the near IR, the Galileo Probe 
Mass Spectrometer (GPMS) data, and ground-based re- 
mote sensing from IR to radio wavelengths. Among the 
detected species, the measured abundances of C, N, S, P, 
Ar, Kr, and Xe (see Fig. [3]) are taken to be representa- 
tive for the convective region below the 1-bar level, the 
outer boundary of our models, as they are not subject 
to non-equilibrium chemistry, cloud formation, or other 
processes that could significantly affect their abundance 
at such low altitudes (jAtreva et al.ll2003l ). Therefore, the 
enrichment of these elements is assumed to be represen- 
tative also for those elements that are still subject to such 
processes, such as oxygen. The Juno orbiter, which is on 
route to Jupiter, is designed to measure the O abundance 
below the cloud level. In Fig. [7] we vary the 0:H ratio 
and calculate the corresponding metallicity of Jupiter's 
atmosphere with the help of Eq. (jG)) , 



E,A^.(N.:H) 



m + MHo(Hc:H) + Y.^ (N^:H) 



(6) 



where fii is the atomic weight of species i. By taking 
into account the measured abundances NiiH with their 



Jupiter models with improved ab initio H EOS 



7 



error bars (mean, minimum, and maximum values) and 
including the non-detected elements {Mg, Al, Ca}, be- 
cause they are abundant in the solar system, with enrich- 
ment factors of 0, 1, and 3x solar, we aim to cover the 
uncertainty in Jupiter's metallicity for a given 0:H ratio. 
The unknown abundances of Mg, Al, and Ca contribute 
an uncertainty to Jupiter's Z\ of ±0.2 x Zq^ and the ob- 
servational error bars of the measured species contribute 
an additional ±0.3 x Zq. 

Assuming Jupiter accreted its volatiles as a result of 
infall of icy planetesimals from cold, outer regions of 
the protosolar cloud which r esembled the interstella r 
medium as we observe it today (jOwen fc Bar-Nunl[l995l) . 
we woul d expect 0:H to be similarly enriched in Jupiter 
as C:H ()Owen et al.l [19991) . i.e. 0:H > 3x solar. For 
this ratio we derive from Fig. [7] a minimum atmospheric 
metallicity of 2 Zq that Jupiter models should satisfy. A 
^ 2.5 X solar metallicity (this work) is consistent with 
an 0:H of 2-4 x solar. A measured 0:H ratio greater 
than 4.5 X solar would not be consistent with the Jupiter 
models presented here. According to the uncertainties 
discussed above, the measured 0:H ratio at the 19 bar 
level sets a lower limit to Jupiter's envelope metallicity 
of 0.8-1.9 Zq (based on Fig.[tl). 

The metallicity of the envelope is directly tied to the 
density of the H-He mixture for fixed (T, P). Our models 
are always warmer (Tc ^ 20, 000 K) than the M08 mod- 
els (Tc ^ 14, 000 K) implying a less dense H-Hc mixture 
that allows to add more metals compared to the MOB 
Jupiter adiabat. This is surprising as non-linear H-Hc 
mixing effects, which have been shown to decrea se the 
density of the mixed phase (jVorberger et al.ll2007[) com- 
pared to the linear mixture that we apply, are already 
included in the H-He EOS of the M08 Jupiter models. 
Non-linear H-He mixing effects remain an appealing pos- 
sibility (see N08) to enhance Z\ in our three-layer Jupiter 
model calculations further. 

We favor the three-layer models over the two-layer 
models as they allow to reproduce J4 without invoking 
the hypothesis of deep-zonal winds because the penetra- 
tion depth of such winds might have to be much deeper 
(< 0.96 i?j) than physic ally allowed in the presence of 
convective flow motions ()Liu et al.ll2008() . If deep-zonal 
winds indeed exist and require a correction of J4 as pro- 
posed in MOB, the atmospheric metallicity of our Jupiter 
models would then rise up to 4.5 Z© and the ice-rock core 
mass up to 12 M©. 

The three-layer model framework is am e nable to ad- 
ditional variations. iFortnev fc NettelmannI (|20100 modi- 
fied the composition of the core by adding envelope ma- 
terial to the core region mimicking a diluted core. For 
extremely diluted cores with rock mass fraction below 
20% in the central region, they found a 50% enhance- 
ment of the resulting atmospheric metallicity of Jupiter. 
Thus the relatively low outer envelope metallicities of our 
models suggest that diluted cores, possibly from core ero- 
sion, may be a better assumption for Jupiter than pure 
rock cores. 

5. CONCLUSIONS 

We have aimed to push three-layer Jupiter models with 
our improved H EOS (H-RE0S.2) in the direction of 
largest possible outer envelope metallicity Z\. We find 
•^1 % 3.0 Z0, corresponding to an 0:H enrichment fac- 



tor < 5. However, we prefer our models with Z\ — 2.0- 
2.5Zo, corresponding to 0:H=l-4 OiHgoiar, because they 
do not require to push all constraints to their \a uncer- 
tainty limits. This increase in the envelope metal con- 
tent compared to our earlier Jupiter models based on the 
H-REOS.l hydrogen EOS arises from about 2% higher 
pressures in the 0.2—1 gcm~^ and 5,000-10,000 K re- 
gion of the H EOS, where the ab initio simulations are 
challenging since dissociation and ionization occurs. The 
resulting Z\ values in our Jupiter models depend on the 
depth of the assumed transition to a helium-rich inner 
envelope. The imposed constraint Z\ > 1 Zq requires 
a transition at Pi -2 > 4 Mbar. This threshold rises to 
Pi -2 > 7 Mbar if fourth order terms in the rotational 
perturbation of Jupiter's gravity field and shape are in- 
cluded. Future investigations of the internal structure 
of Jupiter, and also of Saturn, should include a conver- 
gence check of the models with respect to the treatment 
of rotation. 

Our two-layer mo del calculated with ou r EOS and our 
code differs from the lMilitzer et al.l (|200Bf ) Jupiter model. 
We find a larger envelope metallicity (4.5 Zq) that en- 
ables us to switch over to a three-layer model where the 
envelope is divided into a metal-poor outer and a metal- 
rich inner part to the expense of the core mass. In a 
low metallicity envelope (the MOB Jupiter model), this 
degree of freedom does not exist. The origin of the dif- 
ferences between these two-layer Jupiter models remains 
unexplained. 

For our three-layer models we calculate a cooling time 
of 4.4-5.0 Gyrs, an asymptotic frequency spacing vq = 
156 ± 0.03 /^Hz for global oscillation modes, and a nor- 
malized moment of inertia A = 0.276 ± 0.04%. These 
models are consistent with the available observational 
constraints. The Juno mission will be very helpful for 
further constraining Jupiter's interior by measuring the 
deep atmospheric oxygen abundance, higher-order grav- 
itational moments, and the moment of inertia. 
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APPENDIX 



METHOD OF CALCULATING THE ENTROPY 

The functions P{p, T) and u{p, T) contain the full thermodynamic information. We first re-construct the free energy 
F and then use the Gibbs-Duhem relation to convert to entropy. It allows to derive the specific entropy s(p, T) with 
an offset sq with respect to a reference state (po,To). The offset is unknown but constant for a given reference state. 
First, F is derived with an unspecified offset Fq = F{po,To) by integration of the total differential d{F/T) along some 
path in p — T space. F is related to the entropy by definition, F = U — TS, where U = Mu and S = Ms are the 
extensive internal energy and the extensive entropy, respectively, and AI = Vp is the mass contained in volume V. 
With the reference state (Tp, Vb) and the unknown entropy offset so we can write 

The term in parenthesis in Eq. (jAip is the solution to the line integral 

rT,V 



F{T',V') 
T' 
To.Vo 



(A2) 



and independent of the chosen path of integration if the EOS data are thermodynamically consistent. We choose lines 
of constant density and of constant temperature as our path of integration. While this is a choice of convenience it 
has the drawback of requiring a large scale EOS in T — p space off the planetary adiabat. With dU = —PdV + TdS 
we have 

dT ' dV T ^ ' 

After switching back to [p, r)-space we can write 

JTo,po Jpo P ^0 J To ^ 
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After integrating the right hand side of Eq. (|A4[) we can numerically calculate the entropy 

u{T, p) 1 rT.p 



T 



M 



To, pa 



,F(T'.p') , 

d ^ j,r + So . 



(A5) 



If the underlying EOS is a mixture of different components, the entropy derived from Eq. ()A5|) implicitly includes the 
mixing terms, for instance the ideal entropy of mixing if the EOS is an ideal mixture. The constant of integration sq 
may depend on composition, but not on the initial condition of the planetary adiabat. If a mixture is kept constant 
during a planet's evolution, the offsets sq of planetary adiabats through different outer boundaries cancel each other, 
so that the real entropy difference between two adiabats is then known. 

MOMENT OF INERTIA 

The moment of inertia / of a rigid body with spin-angular velocity w and axis of rotation c„ through the body's 
center of mass is 

/ = ^ p{r) (c„ X r)^ , (Bl) 

where the integral is taken over the volume V of the body. If we choose the orientation of the body such that runs 
along the z-axis, then (e,^ x r)^ = siv? § in spherical coordinates. Because for an isolated fluid planet the assumption 
of hydrostatic equilibrium implies symmetry around the axis of rotation, Eq. (jBip becomes 



I = 2tt 



^0 



dMr pir, d) sin^ 



(B2) 



where R is the planet's radius. If the planet rotates rapidly such as Jupiter, R depends on the latitude i9; in particular, 
the equatorial radius i?oq = Rii9 = 7r/2) is larger than the polar radius R{-d = 0). Therefore, in order to calculate 
the two-dimensional integral the figure of the planet must be known. The figure is the shape at an iso-bar surface, 
equivalently to an equipotential surface. For the solar system gi ant planets, this surface is taken to be the 1 bar 
pressure level by convention. We have used the Theory of figures (jZharkov fc Trubitsvnl[T978| ) to determine Jupiter's 
figure in terms of the figure functions s„(Z) which give the shape of equipotential surfaces 



riil,l9) ^ I i I + J2 S2nil) P2n{cOsd) 



(B3) 



If the new radial coordinate I is fixed, the total potential remains unchanged under variation of latitude The 
functions P2n (cosi?) are the Legendre polynomials. Only even expansion coefficients occur because of the assumed 
symmetry between northern and southern hemisphere. The sum in Eq. (jB3p is truncated after the third order term 
{n = 3) because that is sufficient to calculate the gravitational mo ment Jg. We abbreviate the right hand side of the 
definition ri{l,'d) by r; — l{l + ^). By inserting Eq. ()B3p into Eq. (|B2|) and changing the independent variable from 
r to I using 



dr/dl = 



1 + 



and 1{P — 1 bar) — i?ni, the mean radius, we can finally compute /, 



I = 2tt dd sin^ d / dl 

/o "'0 



In case of spherical symmetry, ~ ^^'^ ds2n/dl = 0. 



dl 



(B4) 



(B5) 
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Figure 1. Hydrogen isotherms of difi'crent EOS; solid: H-REOS.2, dashed: H-REOS.l, and dotted: SCvH-i EOS. Inset: Tlie relative 
deviation of H-REOS.2 with respect to H-REOS.l for selected isotherms (5000 K: solid, 8000 K: dotted, 10000 K: dashed.) 
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Figure 2. Relative difference in density (solid) and temperature (dashed) between the Jupiter adiabats ealeulated with H-REOS.l and 
H-REOS.2, for a H-He mixture with Y = 0.275 (Jupiter average). The improved simulation of hydrogen dissoeiation causes systematically 
lower densities in the outer ~ 16 of Jupiter. Numbers denote the mass in Earth masses between the pressure levels (in GPa) of 0.1-10, 
10-100, 10-1000 as indicated by the arrows, and further up to 40 Mbar close to the core mantle boundary. 
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Figure 3. Jupiter interior models for different transition pressures Pi_2 at the layer boundary. Top: core mass; middle: outer envelope 
mctallicity scaled to the solar value Zgaia,i = 0.015; bottom: inner envelope metallicity. Curves are with the improved H EOS (H-REOS.2) 
and open symbols are models calculated with our former H-REOS.l. Among these, solid curves and triangles are for J4/10~* = —5.89, 
while dashed curves and diamonds are for 74/10-" = -5.87. The thick (thin) curves are calculated with the Theory of Figures to 3rd (4th) 
order. Arrows indicate t he sh ift of the models Jll-4a and Jll-8a [filled circles) if the observational Icr error bars of Y, Yi, Ti, and J4 are 
applied as described in S I3.1I Filled diamonds are measured atmospheric abundances in solar units placed arbitrarily on the x-axis. The 
dotted line is a guide to the eye for a (minimum) outer envelope metallicity of 2x solar. The new models with H-REOS.2 are consistent 
with this constraint, but H-REOS.l based models were not. The dot-dot-dash line in the bottom panel shows the scaled mass coordinate 
mi_2 := m{Pi—2), which is within the line thickness independent on the model assumptions. It can be used to calculate the mass of heavy 
elements in the envelopes Mz^env = Zi X (Mj — mi-2) + Z2 x (mi_2 — M^) . 
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Figure 4. Two-layer models on their way to meet Jupiter's observed J2 at 146.97 X 10~^ (vertical dotted line), starting at low envelope 
metallicity and moving in the direction as indicated by the arrows. The converged solution of our calculations is obtained through variation 
of Z. It has Z = 4.33 Zq and a core mass of 21.1 Afg) if the core is mostly ice (dashed), or Mc = 16.5 Mgj if made of rocks (solid), 
respectively. Square with error bars: the M08 Jupiter model. 




Figure 5. Effective temperature (solid) and mean radius {dashed) of tfie homogeneously cooling Jupiter (interior model Jll-4a). The 
square indicates the observed T^fj at present time, whose observational error bar of 0.3 K is not resolved in this figure. 
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800 600 400 200 

Ti (K) 

Figure 6. Jupiter's evolution including angular momentum conservation, the corresponding change in the energy of rotation, and solar 
irradiation that increases with time. Upper panel: period of rotation; middle panel: normalized moment of inertia A for the surface 
temperatures Ti = 170-800 K that define the internal structure during the evolution. Lower panel: Map of Ti onto time. 
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Figure 7. Mctallicity Zatm in Jupiter's deep atmosphere and outer envelope in dependence on the 0;H ratio in solar units l|Loddersll2003l ). 
Thin solid: using the mean values of the measured species (C,N,S,P,Ne,Ar,Xe,Kr); solid (thick solid): including in addition Ix (3x) solar 
abundances of (Mg,Al,Ca); thick dashed: same as thick solid but using the upper limits of the measu red abundances; thin dashed: same as 
thin solid but using the lower limits. The measured 0:H abundance at 19 bars I IAtreva et al.|[2"003) is indicated (vertical position has no 
meaning) . 
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Table 1 

Jupiter model Jll-4a. 



TO 

(Me) 


P 
(GPa) 


I 

(i?m) 


T 
(K) 


P 

(gcm-3) 


S2 


S4 


S6 


317.833802 


l.OOOOE-04 


1.000000 


170.0 


1.6970E-04 


-4.49641E-02 


1.98582E-03 


-1.01995E-04 


227.158225 
227.155808 


3.9999E+02 
4.0001E+02 


0.735537 
0.735532 


8995.9 
8996.0 


1.3253E+00 
1.4332E+00 


-3.53291E-02 
-3.53290E-02 


1.09438E-03 
1.09437E-03 


-3.82050E-05 
-3.82046E-05 


7.567754 
7.567666 


4.2729E+03 
4.2729E+03 


0.114019 
0.114018 


20041.2 
20041.2 


4.3140E+00 
1.8815E+01 


-1.00806E-02 
-1.00805E-02 


8.05449E-05 
8.05447E-05 


-6.72388E-07 
-6.72385E-07 



Note. — This table is published in its entirety in the electronic edition of the ApJ. A portion is 
shown here for guidance regarding its form and content. The columns show the internal profile of 
the Jupiter model Jll-4a in terms of the mass coordinate m, the pressure P, the radial coordinate I 
scaled by Jupiter's calculated mean radius = 10.95517-Rei the temperature T, the density p, and 
the figure functions S2-S6, see Appendix B. 



